Nonlinear dynamics of beta induced Alfven eigenmode driven by energetic particles 
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Nonlinear saturation of beta induced Alfven eigenmode, driven by slowing down energetic particles 
via transit resonance, is investigated by the nonlinear hybrid magnetohyrodynamic gyro-kinetic code 
(XHMGC). Saturation is characterized by frequency chirping and symmetry breaking between co- and 
counter-passing particles, which can be understood as the the evidence of resonance-detuning. The scal- 
ing of the saturation amplitude with the growth rate is also demonstrated to be consistent with radial 
resonance detuning due to the radial non-uniformity and mode structure. 



Apart from Toroidal Alfven Eigenmodes ^ (TAB), the 
existence of Beta induced Alfven Eigenmodes [J] (BAE) 
in the kinetic thermal ion (KTI) (j, |j] gap of the shear 
Alfven wave (SAW) continuous spectrum is also widely 
recognized to be a concern for the good confinement of en- 
ergetic particles (EP) in fusion plasmas [2, 5]. BAEs are 
particularly important in the study of low frequency fluc- 
tuations of the shear Alfven wave (SAW) spectrum, since 
they can be excited by both fast ions (at long wavelengths) 
as well as by thermal ions (at short wavelengths) la, iTfl- 
In addition to Alfven eigenmodes (AE), which are normal 
modes of the thermal plasma, strongly driven energetic par- 
ticle continuum modes [8] (EPM) may be also excited at 
the characteristic frequencies of the EPs, in the presence of 
a sufficiently intense EP free energy source. 

The different nature of AEs and EPMs suggests the ex- 
istence of two regimes in the nonlinear dynamic evolution 
of a single toroidal mode number (n) coherent SAW driven 
by EPs. Near marginal stability, when the system is weakly 
driven and EPMs are not excited, the nonlinear dynamics 
is the same as in a uniform system and wave saturation 
can occur when wave -particle trapping flattens the parti- 
cle distribution in the resonance region [9]. Meanwhile, 
other physical mechanisms, such as Compton scattering off 
the thermal ions llioll and mode-mode couplings, enhanc- 
ing the interaction with the SAW continuous spectrum by 
nonlinear frequency shift [ Uli \12A . can also be important, 
depending on the parameter regimes. All these phenom- 
ena are local and independent of the radial mode struc- 
ture. However, when the system is strongly driven and 
EPMs [Sl] are excited at the EP characteristic frequencies, 
strong EP transport occurs in avalanches [13]. This phe- 
nomenology is strictly related with the resonant character 
of the modes, which tend to be radially localized where the 
drive is strongest, and with global readjustments in the EP 
radial profiles. In the strong nonlinear regime, local and 
global effects must be treated at the same footing in order 
to capture the crucial role of radial non-uniformity. 

In this letter, we investigate the transition regime in be- 
tween the two limiting conditions described above, in or- 
der to illustrate the effect of equilibrium non-uniformity as 
the mode drive is increased above marginal stability. In 



particular, we analyze AEs resonant excitation below the 
EPM threshold and choose BAE as specific case because: 
(1) rich phenomena are observed in experiments [2, 5]; (2) 
a single poloidal harmonic is dominant; (3) the mode is 
very localized, so that finite mode width effects and finite 
interaction length are easily illustrated. Here, we focus on 
EP nonlinear dynamic behaviors and their consequences on 
BAE saturation. Thermal ion nonlinear dynamics may also 
be important, but it will be analyzed in a separate paper. 

In our work, BAE driven by EPs are studied with 
the extended version of nonlinear magnetohydrodynamic 
(MHD)-Gyrokinetic code HMGC LMl (XHMGC). These 
extensions include both thermal ion compressibility and 
diamagnetic effects [15], in order to account for ther- 
mal ion collisionless response to low-frequency Alfvenic 
modes driven by EPs (e.g., kinetic BAE (KBAE) 11611 ) and 
finite parallel electric field due to parallel thermal electron 
pressure gradient. Previous linear simulation studies show 
that BAE/KBAE can be destabilized by EPs via resonant 
wave particle interactions [|l6l llVh . Here, we ignore dia- 
magnetic effects of kinetic thermal ions by assuming uni- 
form thermal ion density and temperature profiles, but keep 
kinetic thermal ion compression effects, in order to cor- 
rectly describe the KTI gap. Meanwhile, we consider non- 
linear wave-particle interactions with EPs, but neglect both 
mode-mode coupling effects and nonlinear kinetic ther- 
mal ion response. The purpose is to isolate the nonlinear 
physics due to EPs, which do play dominant roles for suf- 
ficiently strong drive. 

We investigate a BAE mode localized around the ratio- 
nal surface at (r/a) ~ 0.5 with safety factor q = 2, 
dominated by toroidal mode number n = 2 and poloidal 
mode number m = 4 harmonic. The BAE-SAW contin- 
uum accumulation point frequency can be calculated from 
kinetic theory ^ and is given by cOcap = 0.127/t^, where 
ta = Rq/Va is the Alfven time, with Va the Alfven 
speed on axis and Rq the major radius. The mode is 
driven by EPs with an isotropic slowing-down distribu- 
tion and characteristic transit frequency at the birth energy 
iOt^a. = {2Eo/mHy^V{qRQ) ^ 0.2/r^. The mode is 
excited at loq ~ 0.114/r^, within the KTI gap, and the 
linear growth rate is ^^ — 0.006/r^. As the mode grows 
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FIG. 1. Upper panel shows the amplitude evolution (black line) 
and frequency evolution of the mode (blue line), given by the 
maximum in the wavelet transform spectrum. Lower panel shows 
the wavelet transform spectrum for the frequency evolution. 



and saturates, the mode frequency is chirping up quickly, 
as shown in Fig. [H because of EP radial redistribution and 
the corresponding nonlinear change in the BAE dielectric 
response. To show this, we use the general fishbone like 
dispersion relation, iA{u}) = 6Wf + 6Wk 10, H], with 

the generalized inertia term iA ~ —\/i-^cap ~ ^"^ I^a for 

BAE yfl, oja = VaUqRo) and va being the Alfven fre- 
quency and speed, respectively, and 5Wf and 5Wk repre- 
senting fluid and kinetic contributions to the potential en- 
ergy [4, 6]. This allows us to write the nonlinear frequency 
shift as 
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having separated linear (L) and nonlinear (NL) contribu- 
tions to the potential energy. Equation ([T]l shows that up- 
ward frequency chirping is expected, due to EP radial re- 
distributions {M.e5WkNL > 0, ImSWkNL < fiH). 
This result stifl holds when iMeJVFfcATi | < \lm6WkNL\, 
since - in the present case of slowing down EPs - 
d^„Im5WkL > 0. 

Due to the preferred direction in mode frequency chirp- 
ing, the symmetry between co-passing (v\\ > 0) and 
counter-passing (v\\ < 0) particles is also broken, as visi- 
ble in Fig. 121 showing the wave-particle power exchange 
in both linear (upper panel) and saturation phase (lower 
panel). Results show that the interaction is predominantly 
driven by transit resonance. Both co-passing and counter- 
passing particles are giving energies to the wave in the lin- 
ear phase, with almost symmetric response in vu space. 




FIG. 2. Wave-particle power exchange in the EP (/i, wii) plane 
during linear phase (upper panel) and saturation phase (lower 
panel), where /i is the magnetic moment (vu and /i are normal- 
ized w.r.t. \/EQ/niH and Eq/^Iq with Qq — ei3o/(m//c), re- 
spectively). The positive sign in the color bar corresponds to 
ions pumping energy into the wave. Solid and dashed lines are 
passing-trapped boundaries, calculated for the inner and outer 
limits of the mode radial width. 



Different fine structures in the resonant region are due to 
different orbits of co- and counter-passing particles. How- 
ever, in the saturation phase, co-passing particles are still 
driving the mode, while counter-passing particles do not. 
This is evidence of resonance detuning, explained as fol- 
lows. The passing particle resonance condition is 



oj = iujt + {nq — m)aujt 



(2) 



where a = sgn{v\\), q is the orbit average of q, q = 
{2tt)~^ <f qdO and the transit frequency is defined as oot = 
2ii/ § dO /9. Here we have used Clebsch toroidal flux 
coordinates (ipjO,^), with the equilibrium magnetic field 
B = Vtp X V^, "0 the magnetic flux and the poloidal an- 
gle 6 describing the position along B. The resonance con- 
dition, Eq. ([21, can be expressed as = 0; i.e., stationarity 
in the wave-particle phase. For the transit resonance under 
investigation, we can define the resonance detuning as 



A0 ~ (nqrcr H 7rdrUJt)(^t'^r 



^oj 



(3) 



where (jr = dq/dr and ujq = Icot + {nq — m)(JLOt. In 
the present condition, the resonance-detuning, due to ra- 
dial particle displacement, is dominated by the first term in 
parenthesis. Thus, for ci; > (upward frequency chirping), 
particles that are transported out, while transferring energy 
to the wave, more easily maintain the resonance condition 
for v\\ > (co-passing) than v\\ < (counter-passing). 
This explains the results of Fig. |2] as well as of Fig. [3] 
where the relative magnetic flux averaged EP fluctuating 
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FIG. 3. Radial dependence of the magnetic flux averaged EP fluc- 
tuating distribution, integrated over /i and on a narrow interval 
Auii near resonance. Upper panel shows the result for co-passing 
particles. Lower panel shows the result for counter-passing parti- 
cles. The radial mode structure with n = 2 and m = 4 is shown 
in the background for reference. 



distribution function, 5n/nQ, is shown vs (r/a), after in- 
tegration in fj, and over a narrow interval Aw|| near reso- 
nance. From these results, it is clear that the wave interac- 
tion with co-passing particles is lasting longer and cover- 
ing a broader radial domain, of the order of the radial mode 
width. 

Equation (jSj also provides useful insights for the analy- 
sis of phase-space structures shown in Fig. |4l There, test 
particle motion in the phase-space is illustrated for both 
CO- and counter-passing particles. For v\\ > 0, the wave- 
particle phase Q is essentially constant, i.e. the G ~ 
resonance condition is effectively maintained through sat- 
uration. This occurs because of Eq. ([T), showing that 
frequency chirps up due to the effect of EP radial trans- 
port. Frequency chirping, in turn, is connected with phase- 
locking (0 ~ 0), for this is the condition that allows parti- 
cles to be most efficiently transported outward, while driv- 
ing the mode. This condition is similar to that, which un- 
derlies EP transports in avalanches [13], due to EPM f^. 
In the present case, however, the finite BAE radial struc- 
ture determines the finite interaction length and mode sat- 
uration is reached when particles are pumped out of the ef- 
fective mode radial domain. For counter-passing particles, 
meanwhile, Fig.|4]shows rapid wave-particle phase mixing, 
as expected from Eq. (|3]l, and the corresponding formation 
of phase-space island structures near resonance, which ex- 
plain the reduced drive for v\\ < in Fig.l2l lll8L[l9ll . This 
result, which crucially depends on radial non-uniformity 



of the system and applies for 
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Wo 



> 



7i, is different 



from the well-known case of nonlinear saturation due to 
wave -particle trapping [l20l - l23n . It is also different from the 
adiabatic frequency chirping of hole-clump structures in 
phase-space II24II . since, in the present case, \uj\ > 'J'lScJ- 
Eq. (|3]l) and the frequency sweeping is non-adiabatic | 
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FIG. 4. Phase space structures in (6, — P^) (6 and P^ are nor- 
malized w.r.t. 2iT and a\/mEo with a the minor radius) phase- 
space for co-passing particles (upper panel) and counter-passing 
particles (lower panel) at saturation, where & is sampled at6 = 0. 



This is a crucial element for maximizing EP transport, for a 
broader region of the phase-space can be affected by radial 
redistribution, as shown in Fig. |3]and|4]by the t;|| symme- 
try breaking. Meanwhile, considering Eq. (H), this explains 
why frequency chirping is connected with the nonlinear 
evolution of the mode amplitude, as shown in Fig. [T] and 
often observed in experiments [25]. 

This interpretation of simulation results is further sup- 
ported by the scaling of the saturation amplitude with the 
growth rate, which is shown in Fig. |5] and can be under- 
stood theoretically by analyzing the nonlinear pendulum 
about the stable fixed point by radially localized perturba- 
tion. The corresponding wave -particle resonance is then 
described by 



dY/dr 
dX/dT 



X 



-AoY/{l + X' 



(4) 



Here, X = \kr\x/{nq') and Y = ylkrl/inq') are the 
rescaled variables for x = nq'Ar, y = @ — ©q, Oq begin 
the fixed point and r = ujtt. We assumed = Aq/{1+X^) 
to investigate the effect of finite characteristic radial mode 
width, i.e. k^^, where Aq is in general proportional to the 
mode saturation amplitude. Equation dU is separable and 
can be trivially integrated by quadratures. As a result, the 
period is given by T = 271 Aq /F, with 

{2/TT){l + Z^X^)dz Y' 



V(l+X„72)-Z2(l + ^2^2/2) 



(5) 
The integral in Eq. (|5]l is a function of Xq, the normal- 
ized maximum radial particle excursion, and is easily com- 
puted by asymptotic expansion, yielding 1 for Xq — )■ and 



s. 



o 




FIG. 5. The saturation amplitude of the fluctuating scalar poten- 
tial vs. the linear growth rate at different values of EP density. 
The "o" are simulation results. The black, blue and red lines ai'e, 
respectively, denoting oc (7/0;)^, (7/0;)'^, (7/0;) for references. 



oc Xq for Xq ^ 1. The characteristic rate for a particle 
to undergo a radial excursion Xq = \kr\Aro, is then given 
by Tatl = A„' ujtT{Xo). If Xq = |A;,.|Aro < 1, as as- 



sumed in the nonlinear saturation of TAE modes H, the 
nonlinear dynamics is the same as in a uniform system and 
wave saturation can occur only when the particle distribu- 
tion function in the resonance region is flattened, so that the 
nonlinear drive is significantly reduced and brought back 
to threshold. This occurs for ^nl/^i ~ Ih/^t ~ ^n , 
yielding the well known estimate Aq cx {'-^l/ojiY 120- 
2211 ■ When Xq = |A;r|Aro > 1, saturation occurs be- 



cause particle get off resonance by radial detuning, i.e., af- 
ter a radial displacement of the order of the mode width, 
as noted above and observed in Il26ll for the first time. 
Similar to the uniform case, saturation is expected when 

-iNL/oot = AI'^T{Xq) a AI'^/Xq « -iL/uJt. Mean- 
while, noting that Y ^ \kr\/{nq'), Eqs. (|4]l also give 

1 /9 

Xq oc {Aq /XQ)\kr\/{nq'). Thus, the saturation con- 
dition becomes Aq oc {'^L/LoY{krY /{nq'Y. By further 
increasing the drive, the radial scale l/{nq') of the BAE 
mode will become increasingly more important, till the 
saturation amplitude is independent of (7^/0;). For BAE 
modes, the transition between these two regimes is ex- 
pected to occur for Xq ~ 1, i.e. for Aq ~ {nq'Y /k^, 
where the radial wave number estimate for the BAE short 
scale radial structure is \kr\/{nq') oc (Jw/u;)"^^^ [3], i.e. 
\kr\/{nq') oc {-iL/^)~^''^ for 7l > \M.e5bj\. The ex- 
pected scaling is then given by Aq oc (^l/uj)^. Figure |5] 
illustrates the transition between these various regimes as 
7l/li; is increased, further demonstrating the important 
role of plasma non-uniformity and radial mode structures 
for mode saturation and EP transport processes. 

In summary, we have found that nonlinear saturation of 
BAEs, driven by slowing down EPs via transit resonance, 
is characterized by upward frequency chirping and v\\ sym- 
metry breaking between co- and counter-passing particles, 
which can be understood as the the evidence of resonance- 



detuning. Upward frequency chirping is the preferred con- 
dition for maximizing EP radial transport while they drive 
the mode, and saturation eventually occurs when the par- 
ticle radial displacement is of the order of the radial mode 
width. The scaling of the saturation amplitude with the 
growth rate further demonstrates the role of radial reso- 
nance detuning, due to the radial non-uniformity and mode 
structure. Both EP phase space behaviors as well as nonlin- 
ear dynamics and mode structures demonstrate the crucial 
roles of non-uniformities and geometry of the system. 
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